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Abstract 

Experiments with trapped atomic gases have opened novel possibilities for studying 
the evolution of nonequilibrium finite quantum systems, which revived the necessity 
of reconsidering and developing the theory of such processes. This review analyzes 
the basic approaches to describing the phenomena of equilibration, thermalization, 
and decoherence in finite quantum systems. Isolated, nonisolated, and quasi-isolated 
quantum systems are considered. The relations between equilibration, decoherence, 
and the existence of time arrow are emphasized. The possibility for the occurrence of 
rare events, preventing complete equilibration, are mentioned. 
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1 Trapped cold atoms 



1.1 Introduction 

In recent years, there has been a remarkable progress in experimental studies of cold trapped 
atoms, which has been accompanied by intensive theoretical investigations. Numerous liter- 
ature can be found in the books [1-4] and review articles [5-19]. Experiments with trapped 
atoms possess several important features making such experiments one of the most conve- 
nient tools for studying the properties of many-body systems. First, trapped atomic systems 
allow for a high degree of tunability, which makes it possible to vary their parameters in a 
rather wide range. Second, the variation of system parameters can be realized very quickly, 
thus, putting a system in a nonequilibrium initial state and studying the system evolution. 
Third, the diluteness of the cold atomic gases and exceptionally low temperatures result in 
rather long timescales of dynamical effects (typically of the order of milliseconds or longer) 
allowing for a very precise time resolution of nonequilibrium processes. Finally, trapped 
atomic systems can be well isolated from environment. Therefore, it is possible to study the 
quantum dynamics of finite systems under different levels of isolation and to analyze the 
influence of the latter on the system equilibration, decoherence, and thermalization. 

This review is devoted to exactly the latter problems. In the present section, the most 
interesting experiments and computer simulations, considering the dynamic equilibration 
effects of nonequilibrium trapped atoms, are surveyed. While the following sections are 
devoted to the general theoretical description of such processes. In order to avoid ambiguities, 
the review tries to use accurate mathematical formulations. 

Throughout the paper, the system of units is employed, where the Planck and Bolztmann 
constants are set to unity (h — 1, ks = 1). 

1.2 Experiments with trapped atoms 

In the experiment of Greiner et al. [20] 87 Rb atoms were loaded into a three-dimensional 
optical lattice. Varying the ratio U / J of the one-site interactions U to the hopping parameter 
J, the transition between the superfluid Bose-condensed phase and the incoherent Mott 
insulator phase was studied. When atoms were initially in the Mott insulating phase with 
the ratio U/J > 5.8z , with z being the number of nearest neighbors, and then this ratio 
was quickly reduced to U/J < 5.8^ , corresponding to the superfluid phase, the coherence 
was restored during the timescale of 1/ J . In the opposite process of quench from the Bose- 
condensed superfluid to the Mott insulator state, the equilibration was accompanied by 
collapses and revivals of the coherent fraction on the shorter time scale of 1/U [21]. The 
equilibration, after crossing the superfluid-insulator phase transition line, evolved through 
the inhomogeneous state with alternating regions of incoherent Mott insulator phases and 
coherent superfluid phases [20]. 

The dynamics of pair condensate formation in a strongly interacting Fermi gas of 6 Li, 
close to Feshbach resonance, was studied by Zwierlein et al. [22]. The relaxation time, after 
a rapid magnetic field ramp across the Feshbach resonance, was of the order of 100/ E F . Here 
Ep ps w (3iV) 1//3 is the gas Fermi surface energy, wo is the geometric mean of the trapping 
frequencies, and N is the total number of atoms in the trap. After a rapid ramp, molecular 
Bose condensate co-existed with unpaired Fermi atoms. 



4 



Trapped quasi-one-dimensional arrays of Bose gases, each containing from 40 to 250 
87 Rb atoms, were prepared in an out-of-equilibrium state by Kinoshita et al. [23]. They 
observed that the system of two initially separated Bose-Einstein condensates oscillates, as 
in a Newton cradle, but does not noticeably equilibrate even after thousands of collisions. 

Spontaneous symmetry breaking in a quenched spinor Bose-Einstein condensate was 
studied by Sadler et al. in a quasi-two-dimensional trap [24]. Nearly pure Bose conden- 
sate of spinor 87 Rb was prepared in the unmagnetized m z = state at a high quadratic 
Zeeman shift. By rapidly reducing the magnitude of the applied magnetic field, the gas 
was quenched to low quadratic Zeeman shift, obtaining conditions that favored the ferro- 
magnetic phase. Spontaneous symmetry breaking occurred through a heterogeneous stage, 
where small ferromagnetic domains of variable size and orientation were intermixed with 
unmagnetized regions. This could be regarded as the phase separation of magnetized and 
unmagnetized states in the process of nonequilibrium phase transition. Spin-vortex defects 
were also observed with high confidence in about one-third of all images containing significant 
ferromagnetism. Some images indicated as many as four vortices. All observed vortices were 
singly quantized, with no apparent preferred direction or circulation. The thermalization 
time was an order longer than the local equilibrium time. 

Nonequilibrium coherence dynamics of quasi-one-dimensional Bose gas of Rb was con- 
sidered by Hofferberth et al. [25]. The initial state consisted of two completely separated 
Bose condensates loaded into a double well, where they started interacting and developing 
the common coherent system. The final state was a thermalized Bose-condensed gas with 
a temperature slightly higher than the initial temperature of separated Bose condensates. 
The thermalization time, depending on the system parameters, was between 6 and 9 ms. 

The spontaneous formation of vortices in the process of nonequilibrium Bose-Einstein 
condensation and the relation of this phenomenon to the Kibble-Zurek mechanism [26,27] 
was investigated by Weiler et al. [28]. 

In the experiment of Hung et al. [29] , a 133 Cs quantum Bose gas was trapped in a two- 
dimensional optical lattice. The lattice depth was quickly ramped to induce the superfluid 
to Mott-insulator phase transition. The global equilibration involved mass transport and 
was found to be much slower than the microscopic tunneling time scale 1/ J. 

By using the high spatial resolution in-situ fluorescence imaging, Bakr et al. [30] and 
Sherson et al. [31] demonstrated that the superfluid to Mott insulator phase transition can 
be studied with single-site and single-particle resolved detection. 

1.3 Numerical simulations with models 

Numerical simulations for different models play the role of computer experiments, providing 
essential information on the relaxation processes in nonequilibrium finite quantum systems. 
In the present subsection, some of such simulations are briefly described. The detailed 
analysis of terminology and correct mathematical definitions will be given in the following 
sections. Many of simulations were done for lattice models. Such models are also of interest 
for cold atom physics, since periodic atomic systems are now easily created by means of 
optical lattices. 

Probably, the first computer simulation of nonequilibrium dynamics for a finite classical 
system was done by Fermi et al. [32] who considered a one-dimensional chain of anharmonic 
oscillators and discovered that it does not equilibrate. But that was a classical system. 
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A quantum two-dimensional anharmonic oscillator, corresponding to the Henon-Heiles 
model [33], was studied by Feingold et al. [34]. They observed that the mean location and 
dispersion of the wave packet, first, tend to quasi-equilibrium values, but later display large 
fluctuations. 

Jensen and Shankar [35] studied relaxation of a finite spin-half quantum chain in a mag- 
netic field in two situations, when the model was integrable or not. They attempted to find 
out if there would be difference in the relaxation of the integrable and nonintegrable cases. 
It turned out that both integrable and nonintegrable quantum systems, with as few as seven 
degrees of freedom, could exhibit thermalization in finite time to the states well described 
by microcanonical ensemble, provided that the initial state was not an eigenstate of any 
considered quantities. Of course, the thermalized state was not an equilibrium state in the 
strict sense, but rather a quasi-equilibrium state, since the overall behavior of a finite quan- 
tum system should be quasi-periodic in time. But that was a quasi-equilibrium state with 
a long lifetime corresponding to recurrence time. An important observation was that the 
equilibration of integrable as well as nonintegrable systems was very similar to each other. 

Rigol et al. [36] and Cassidy et al. [37] simulated the relaxation of one-dimensional lattice 
hard-core bosons. This is an integrable model that equilibrates to a representative Gibbs 
ensemble [38,39], differing from the standard ensembles, similarly to the integrable Luttinger 
model, as has been shown by Casalilla [40] . A rigorous definition of the representative statis- 
tical ensembles will be given in the following sections. Generally, the class of representative 
ensembles includes the standard microcanonical, canonical, and grand canonical ensembles. 
To distinguish between them, the representative ensembles that are different from the latter 
can be called nonstandard. 

Comparing the behavior of integrable and nonintegrable systems for one-dimensional 
[41-43] and two-dimensional [44] lattices of hard-core bosons, Rigol et al. [41-44] concluded 
that integrable systems relax to a nonstandard representative Gibbs ensemble, while the 
nonintegrable systems relax to a standard microcanonical ensemble. Similar behavior occurs 
for one-dimensional lattices of fermions [45]. 

Kollah et al. [46] investigated one- and two-dimensional optical lattices described by 
the boson Hubbard model. They studied the time evolution following a quench from the 
superfluid to the Mott insulator phase. The quench was realized by a sudden change of the 
on-site interaction U under a fixed hopping parameter J. The superfluid-insulator line, in 
the case of unit filling, is characterized by the critical value u c of the ratio u = U/(Jz ), 
where zq is the number of nearest neighbors. For one-dimensional latices, u c = 1.8 and 
for two-dimensional lattices, u c = 4.2 (see review [18]). The quench from the superfluid to 
insulator phase corresponds to the change of u from u < u c to the values u > u c . For large 
values of the final interaction strength u ^> u c , the system reaches a quasisteady state that 
is very different from the standard thermal equilibrium and retains memory of the initial 
state. As the final u decreases, the nature of the steady state changes, and in the region 
where U is comparable to J, that is, u ~ l/^o, the steady state is well approximated by 
a standard thermally equilibrium Gibbs ensemble. After the quench, the system exhibits 
a number of oscillations with a period 2tc/U, which relax to a quasisteady state on a time 
scale 1/ J. Similar results for a two-dimensional lattice of bosons were found by Natu et al. 
[47]. 

The existence of two distinct equilibration regimes seems surprising, given the noninte- 
grability of the boson Hubbard model for any finite values of U and J. This tells us that the 
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equilibration regime and the quasisteady state resulting in this process depend mainly not 
on the integrability of the system but rather on the system parameters. Both, integrable as 
well as nonintegrable quantum systems, depending on their parameters, can relax either to 
a standard statistical ensemble or to a nonstandard one, any of such ensembles being just 
particular cases of the class of representative ensembles. 

The results, similar to the boson Hubbard model were found for a one-dimensional 
fermion Hubbard model by Manmana et al. [48]. The ground-state phase diagram of this 
model, at half filling, is well known. For weak on-site interactions U, lower than the critical 
value, given by U c / J = 2, that is, u c = 1, the system consists of a Luttinger liquid, while 
for U > U c , it is a charge-density-wave insulator. By a sudden change of the interaction 
strength it is possible to switch the initial state between metallic and insulating phases. The 
relaxation goes through an oscillatory stage, with the oscillations of period 27r/U, and tends 
to a steady state on a time scale 1/J. The reached quasistationary state is not thermal 
in the usual sense, since it depends on the initial conditions. The time evolutions, starting 
from different initial states, can be distinguished from each other, so that neither relaxation 
to one and the same equilibrium state nor the usual thermalization occur. The resulting 
quaisisteady state could be described by a representative ensemble with fixed expectation 
values of the powers of the Hamiltonian < H n >. The choice of these additional conditions 
can be justified by the fact that, for a closed system, the set of all powers of the Hamiltonian 
constitute an infinite number of integrals of motion. 

The considered one-dimensional fermion Hubbard model [48], with nearest-neighbor in- 
teractions, is integrable. The authors have also studied the effect of adding to the model 
next-nearest-neighbor interactions, which makes it nonintegrable. It has been found that, 
independently of the model integrability or criticality, the relaxation process does not change 
much and the system always relaxes to a nonthermal quasistationary state that can be de- 
scribed by a representative ensemble with additional conditions of fixed < H n >. When 
different initial states have the same energy and are sufficiently close to each other, the 
observable quantities relax to the same value, i.e., the memory of the initial state in that 
case is lost after the relaxation. However, the general energy distribution is defined by a 
representative ensemble retaining the memory of the initial states. 

A chain of two-level systems interacting through the van der Waals forces, with a strength 
U, was considered by Lesanovsky et al. [49]. The transitions between the energy levels 
were due to the Rabi term, with a Rabi frequency Q. For weak interactions U/Vt <C 1, 
no equilibration was noticed. But for sufficiently strong interactions U/Q > 2, the system 
reached a steady state at the time of order 50/f2. The resulting steady state could be well 
characterized by a grand canonical ensemble. Thermalization also occurs for a finite spin 
system with randomly distributed initial states [50]. 

Integrable as well as nonintegrable finite quantum systems, it seems, both can equilibrate. 
Integrability does not seem to play a crucial role on the structure of the quasistationary 
state. This is despite the fact that integrable and nonintegrable quantum systems display 
different level-spacing statistics and differently react to external perturbations [51]. Though 
integrable systems can equilibrate, but their main difference from nonintegrable systems can 
be in much longer equilibration times. This behavior is contrary to integrable classical finite 
systems that do not equilibrate at all. Nonintegrable classical systems can equilibrate if they 
are chaotic [52,53]. 

The influence of the vicinity of critical points on the nonequilibrium dynamics of closed 
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quantum systems was investigated by Polkovnikov et al. [54-56] who emphasized the uni- 
versality of this dynamics in gapless systems near continuous phase transitions. 



2 Thermalization versus equilibration 
2.1 Quantum statistical ensemble 

The term statistical ensemble is constantly employed in studies of many-body systems. It is 
therefore important to recall its correct mathematical definition. 

A quantum system is characterized by a Hilbert space of microstates J 7 . The probability 
operator measure is given on J 7 by a statistical operator p(t) parameterized by the time 
variable t > 0. Hence p(t) is a positive operator normalized to one, 



where the trace is over J 7 . By definition, a quantum statistical ensemble is the pair {J 7 , pit)}. 
2.2 Algebra of local observables 

On the Hilbert space J 7 , one defines self-adjoint operators A that are called the operators of 
local observables, provided their expectation values 



correspond to the measurable quantities. The set of all such operators forms the algebra of 



2.3 System statistical state 

The set of all expectation values for the operators of local observables is the statistical state 



Trp(t) = 1 , 



<i(t)> = Trp(t)A 




local observables A = {^4}- 




(2) 



An important role in dynamical problems is played by the initial statistical state 




(3) 



which is the set of the expectation values for the operators of local observables, 



(A(0)) = Trp(0)A , 



at the initial moment of time t = 0. 
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2.4 Equilibration of quantum systems 

In a review, devoted to equilibration and thermalization, it is crucial to follow precise defini- 
tions of these processes. In literature, one can meet different descriptions of these phenomena. 
Below, the definitions are given that will be used in the present review. 

A statistical system, characterized by a statistical ensemble {J 7 , p(t)}, equilibrates from 
an initial statistical state Ao if and only if for any A G A there exists a limit 

hm(A(t))=Trp*(Ao)A, (4) 

generally, depending on the initial state Ao- 

If the limit ([4j) exists, one says that the statistical system is in an equilibrium state with 
a statistical ensemble {J 7 , p*(Ao)}. 

This is a general definition of equilibration as such. One also considers another form of 
equilibration, the equilibration on average. This will be treated in the following sections. It 
is worth stressing that the above definition of equilibration requires that the limits @ would 
exist for all operators from the algebra of local observables. When such limits exist only for 
some observables, but are absent for others, it is not, strictly speaking, equilibration. 

Also, when a statistical state appears to be time-independent for a finite interval of time, 
but not for the actual limit t — > oo, this is termed quasi- equilibration. 

2.5 Thermalization of quantum systems 

Thermalization is a more restrictive process, when not merely limit @ exists, but the mem- 
ory of initial conditions is, at least partially, lost. 

A statistical system, characterized by a statistical ensemble {J 7 , p{t)}, thermalizes from 
an initial state Ao if and only if there exists a dense set B = {Ao} of initial states, including 
the given initial state Ao, such that for any A G A the limit 

lim(A(f)) =Trp*i (5) 

t— >oo 

does not depend on Ao G B. 

The set B is termed the attraction basin of the stationary state (JSJ). 

When the thermalization occurs from any available initial state, that is, when the attrac- 
tion basin is the whole set of all admissible initial states, this is called the global thermaliza- 
tion. 



3 Evolution of statistical states 
3.1 Representative statistical ensembles 

To correctly describe the system equilibration, it is necessary to be accurate in classifying 
the limiting stationary states. Each statistical system is characterized by a statistical en- 
semble. It is assumed that the statistical ensemble is defined so that to uniquely represent 
the considered statistical system. This implies that the definition of the ensemble must 
take into account all conditions and constraints that correctly represent the system. This 
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principal point was emphasized by Gibbs [38,39], Tolman [57] and ter Haar [58]. The term 
representative ensemble was introduced by Tolman [57]. 

The practical realization for constructing a representative ensemble can be illustrated for 
an equilibrium system, following the basic Gibbs idea. Let a Hilbert space of microstates J 7 be 
given. To define a statistical ensemble {J 7 ,/)}, one needs to specify the operator probability 
measure, that is, the statistical operator p. In constructing the latter, it is necessary to take 
into account those constraints that uniquely define the system. Of course, the statistical 
operator is to be normalized as 

Trp = 1 . (6) 

An important quantity is the internal energy given by the average of the energy operator H, 
called Hamiltonian, 

TrpH = E . (7) 

In addition, there can be other constraints defining the expectation values of some constraint 
operators Ci, with -z = 1,2,..., 

TtpCi = d . (8) 

The statistical operator of an equilibrium state is defined as the minimizer of the infor- 
mation functional 

I[p] = Tip In p + A (Trp - 1) + (TipH - E) + (3 £ \ (irpCi - Ci) , (9) 

i 

in which Ao, Aj, and /3 are Lagrange multipliers. This is equivalent to the maximization of 
the Shannon entropy under conditions ((BJ-dSJ). The result is the statistical operator 

P = \ e-? H , (10) 
with the normalization factor being the partition function 

Z = Tre~ pH , 

and with the grand Hamiltonian 

H = H + J^^Ci- (11) 

i 

As is clear, when the sole constraint is the normalization condition (JS)) , then the statistical 
operator reduces to the microcanonical form 

p = 1 (Z = Trl) . 

If only conditions (JHJ) and ((7j) are accepted, one gets the Gibbs canonical ensemble. If one 
also defines the average number of atoms, one gets the Gibbs grand canonical ensemble. The 
general situation, with any number of additional constraints, is also termed the grand Gibbs 
ensemble [59,60]. 

The general form of the grand ensemble, with the statistical operator ( 1101) . was introduced 
by Gibbs [38,39]. One also calls this the generalized or conditional Gibbs ensemble [61]. But 
one should not forget that this ensemble was advanced by Gibbs. When deciding on a 
particular type of a Gibbs ensemble, one has to choose such that would be representative 
for the considered statistical system. General mathematical properties and applications of 
representative statistical ensembles have been analyzed in Refs. [13, 62-64]. 
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3.2 Nonequilibrium representative ensembles 



When a quantum system equilibrates, it tends to one of the equilibrium representative sta- 
tistical ensembles. A nonequilibrium state is characterized by a nonequilibrium ensemble 
{J 7 , p(t)}. The temporal evolution of the statistical operator is governed by a unitary evo- 
lution operator U(t), such that 

p(t) = U(t)p(0)U + (t) . (12) 

The evolution is generated by the grand Hamiltonian f TTTj) . and the evolution operator sat- 
isfies the Schrodinder equation 

i ± U(t) = HU(t) . (13) 

This means that the grand Hamiltonian is the evolution generator. 

If one assumes that the evolution is generated by the energy Hamiltonian H, then the 
constraint operators Cj in the grand Hamiltonian (11 ip . defining an equilibrium state, are 
to be the integrals of motion, commuting with H. But if the evolution is generated by the 
grand Hamiltonian (fTTp itself, the constraint operators do not need to be the integrals of 
motion [13,62-64]. 

In particular, the constraint operators can fix the required initial conditions for observable 
quantities in the form 

Q = (&) = Trp(0)A . (14) 

These initial conditions then define the solutions to the evolution equations for the studied 
observables. Generally, the steady states, resulting from sudden changes of the Hamiltonian, 
are not thermal, in the sense that they retain information on initial conditions [65,66]. As 
examples of the cases, where the memory of initial conditions is retained in quasi-equilibrium 
states, it is possible to recall the hard-core bosons in lattices [36,37,41-44], Luttinger liquids 
[40], Hubbard-type models of bosons and fermions [45-48,67,68], free bosonic models [69], 
and one-dimensional fermionic systems with Hall states [70]. Other examples are given 
by the quasistationary states that are the solutions to the evolution equations describing 
the relaxation of average spins for different strongly nonequilibrium spin systems (nuclei, 
electrons, molecules, clusters), whose dynamics essentially depends on initial conditions, 
such as the initial spin polarization [71-80]. 

It is possible to conclude that closed quantum nonequilibrium systems equilibrate to 
quasi-equilibrium states characterized by representative statistical ensembles. In the ma- 
jority of cases the resulting representative ensemble retains the memory of initial states, 
regardless of whether a model is integrable or not. 



3.3 Evolution of observable quantities 

What one needs and measures in experiments is the statistical state, that is, the averages 
of operators of observables. Therefore, equilibration has to do with the temporal behavior 
of such averages. An average of an operator A, because of the evolution law ( fT2l) . can be 
written in two ways: 

(A(t)) = Trp(t)A = TrpA(t) , (15) 

in which 

P EE p(0) , A EE i(0) , 
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and the operator time evolution is given by the Heisenberg form 

A(t) = U + (t)A(0)U(t) . (16) 

The expectation value ( !T5|) is real, as a consequence of the operator of an observable being 
self-adjoint. 

The Heisenberg equation of motion for the operator (fTBT) reads as 



i | A(t) = zU + (t) — U(t) + [A(t), H(t)\ , (17) 

where 

H{t) = U + (t)HU(t) . (18) 

The Hamiltonian, being the evolution generator, varies with time only through an explicit 
dependence on time, if any, 

J PjTT 

t -H(t)^U+(t)-U(t). (19) 

If the considered operator does not explicitly depend on time, then its average evolves in 
time according to the law 

4<i( t )) = ([i( t ), H [t )}) (f = ). (20) 

The evolution equation is to be complemented by an initial condition A(0) =< A(0) >. It 
is the memory of such initial conditions that can be retained in the representative ensemble 
characterizing equilibrated finite systems. 

3.4 Equations in matrix representation 

Let {\n >} be an orthonormal basis in the Hilbert space of microstates J 7 . Define the matrix 
elements for the statistical operator 

Pmn (t) = (m\p(t)\n) (21) 
and for the operator of an observable 

A mn (t) ee (m\A(t)\n) . (22) 
The related matrix elements at zero time are 

Pmn = (m\p\n) = p mn (0) , A mn ee (m\A\n) = A mn (0) . (23) 
Then the average (IT5|) can be represented as 

(t) . (24) 

mn mn 

This average, in general, contains nondiagonal terms. Of course, if one chooses as the 
basis the set of the eigenvectors of A, then the sum would contain only the diagonal terms. 
However, other operators, noncommuting with A, will have nondiagonal terms. When one 
discusses equilibration, one should keep in mind not just one operator but the set of operators 
from the algebra of local observables A. Equilibration, by definition, presupposes that it is 
the whole statistical state < A > that equilibrates, but not merely some of its components. 
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3.5 Pure quantum states 

One often considers equilibration of quantum systems prepared in a pure state described by 
a normalized wave function |^(0) >. With a wave function \ip(t) >, the statistical operator 
of a pure state is 

m = \m)(m\- (25) 

It is an idempotent operator, with the property p 2 {t) = p(t). The wave function satisfies the 
Schrodinger equation 

i j t \m) = H\m) ■ (26) 

Employing an expansion 



one gets the equation for the expansion functions 

d 



i — c n (t) = ^H nm c m {t) , (2£ 



where H mn =< m\H\n >, complimented by an initial condition cq = c(0) and the normal- 
ization condition 

n 

The matrix element H mn can be explicitly dependent on time. 
The statistical operator (125]) enjoys the expansion 

P(*)=X}Pmn(*)NH, (29) 
mn 

in which 

Pmn(t) = C m (t)c* n (t) . 

Average fl2l]) takes the form 

(A(t)) =Y J Pmn(t)A nm = (^(t)|i|^(t)> • (30) 
mn 

The pure-state statistical operator ( 1251) is a particular case of the general statistical operator. 
Because of this, it is sufficient to consider the general case. 

4 Isolated quantum systems 
4.1 Absence of absolute equilibrium 

Let us assume that a finite quantum system can be completely isolated from its surrounding. 
Hence the Hamiltonian H does not explicitly depend on time. Then the evolution operator 
is 

U(t) = e-*™ (¥L = o) . (31) 
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It is admissible to take the set of the Hamiltonian eigenvectors, defined by the eigenproblem 



H\n) = E n \n) , (32) 
as the basis in J 7 . Then the matrix element (1211) becomes 

Pmnif) = Pmrfi m ™ > (33) 

with the transition frequency 

ui mn = E m — E n . (34) 

The operator average fl2~4"l) yields 

(A(t)) = J2pmnA nm e- l ^ t . (35) 



And the equation of motion (|20|) reduces to 

d 
It 



I — (A(t)) = J^PnnAnmUmne-™™* ■ (36) 



It is evident that Eq. (I35I) is a quasi-periodic function of time. Hence, it cannot tend to 
a time-independent stationary state. Any given initial state will reproduce itself after the 
Poincare recurrence time [81-84]. The temporal limit (j4]) does not exist. A finite quantum 
system, completely isolated from its surrounding, does not have an absolute equilibrium in 
the sense of limit 

An exception could be the case, when the system would be initially prepared in a pure 
state, with a given Hamiltonian eigenvector \j >, when p mn = 5 m j5 n j. Then the average 
never changes with time: 

(A(t)) = Ajj (p mn = 5 mj 5 nj ) . (37) 

However, the dynamics of the system with this initial state is not structurally stable. If the 
initial condition slightly deviates from the pure state \j > or there is a perturbation of the 
Hamiltonian, so that the initial condition gives 

Pmn ^mj^nj vPmn ; 

then the average 



(A(t)) = A j3 + J2( 6 Pmn)A 

mn 

again becomes quasi-periodic. 



iuj mn t 



4.2 Possibility of equilibrium on average 

Von Neumann [85] suggested that quantum systems could exhibit a kind of ergodic behavior, 
similarly to that of classical systems, when the time average of the expectation value of an 
observable coincides with its ensemble-averaged value. For this purpose, one considers the 
time average 

(A(t)) = lim - [ (A(t)) dt . (38) 
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For an isolated system, this gives 



) , (39) 

ran 

where 

A(w) = lim - / e~ iwt dt , 

r^oo T Jq 

which gives 



AM 



1 , co = 
, w ^ 



If the eigenenergies of the Hamiltonian are not degenerate, then 



A(u mn ) = lim - I e- 1 ^ dt = 6 mn . (40) 



On usually assumes this nondegeneracy condition, which yields 



{A(t)} =J2pnnA nn . (41) 

n 

Hamiltonians, enjoying nondegenerate spectra, are termed typical or generic. 

Von Neumann assumed that, under some conditions, the right-hand side of Eq. f f4T]) 
would correspond to the average over an equilibrium ensemble. In that sense, Eq. (j4"T|) 
would play the role of an ergodic relation, similar to that in classical systems. 

It may happen, however, that, even when this ergodic- type relation is valid, but system 
fluctuations around the averaged value are so large that it is senseless to talk about equilib- 
rium. In order that it would be possible to state that a quantum systems spends the most of 
its time close to the averaged value, having not so many and not so large deviations from it, 
it is necessary to analyze the strength of the system fluctuations. For an operator A, these 
fluctuations are quantified by the dispersion 



a\ = (A(t))* - (A(t)) , (42) 



which is also called the variance. 

Considering the dispersion, one usually assumes that, in addition to the nondegeneracy 
condition ( 1401) . the nonresonance condition 

lim - / e^ 4 *")' dt = 5 mn 5 kl + 5 ml 5 nk - 5 m J nk 5 kl (43) 
r^oo r Jq 

holds true. If so, it immediately follows that 

<y\ = ^1 \pmnA nm \ 2 . (44) 



One also employs the inequality 



y^Pnn < (max p. 



a-l 
nn J t 
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being valid for any real a > 1. 

Another important point is the Chebyshev inequality stating that for any random variable 
x, with average x and variance cr 2 , and any given e > 0, the probability that x deviates from 
x by more than e satisfies the inequality 

P{\x-x\ > e) < (?~y . (45) 

This inequality allows for the estimation of the operator variance (1441) characterizing the 
deviations from the operator average [86-88]. 

For the operator average < A >, assuming that the experimental resolution 5a of 
measuring the observable, associated with the operator A, is sufficiently small, and that 
max n p nn <C 1, Reimann [87,88] obtains 



(A(t)) - (A(t))\ > 5 A ) < . (46) 



This inequality tells us that, under some conditions, the deviations of the observable quantity 
< A > from the time-averaged value (A) could be relatively small. Hence the system would 
spend its major time near the time-averaged value. 

The assumption that max n p nn <C 1 presupposes a uniform distribution over many degrees 
of freedom, which is equivalent to ascribing to the system a high effective temperature. Such 
an assumption loses its validity at low effective temperatures, when it can be that p nn ~ 1, 
for instance, when there is Bose-Einstein condensation. 



4.3 Eigenstate thermalization on average 

If the initial statistical operator corresponds to a fixed pure state \ j > that is a Hamiltonian 
eigenstate then, according to Eq. (|37|) . the observable quantities do not change in time, 

(A(t)) = A 33 = (I) , (47) 

which is evident, since the Hamiltonian eigenstate is a stationary state. It is reasonable 
to expect that there could exist a set of states around a given eiegenstate from which the 
system would equilibrate on average to a state close to that eigenstate. 

Let us fix the eigenstate \j >, whose energy is Ej. And let us consider an energy shell 

Ej = {E n : \E n -Ej\< AEj} (48) 

of energies deviating from Ej not more than by AEj. The energy shall is characterized by 
the indicator function 

If we take for an initial condition an eigenvector with the energy in the energy shell (148|) . 
then we have 

Pnn = Pnj£,j{E n ) , (50) 
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with the normalization condition 



J^Pnn = ^Pnjij{E n ) = 1 . 



n 



The variation of the diagonal matrix elements of an operator, in the range of the energy 
shell, is smaller than 

AAj = max A nn — min A nn , (51) 



which is a non-negative quantity. The values A nn lie in the vicinity of the average 



It is always possible to choose so narrow energy shell that 

AA 



< 1 . (53) 



Then the elements A nn are almost constant inside the chosen energy shell, so that one can 
invoke the theorem of average 

J2PnME n )A nn ~ A^PnMEn) = A* . (54) 

n n 

Therefore, for any normalized weight p n j, we have 



(A(t)) ~ A* . (55) 

Equality ( 155]) implies that there exists an attraction basin in the set of initial states, from 
which there happens an approximate equilibration on average to the value A*. Since, under 
condition (|53|) . the weight p n j can be taken as an arbitrary normalized distribution, it can 
be chosen to be a simple uniform distribution 

This allows us to rewrite average (|52p in the form 

Aj = ^ ^ Pnj£,j(E n )A nn , (57) 

n 

which corresponds to the Gibbs microcanonical ensemble. Such an equilibration can be called 
the eigenstate thermalization on average. Different variants of this type of thermalization 
have been considered in Refs. [89-92]. 

It is worth emphasizing that this thermalization is approximate, it happens on average, 
initial states are assumed to be pure and sufficiently close to a fixed eigenstate, the energy 
shell has to be rather narrow, with the width of this shell depending on the considered 
observable, and also the width can depend on the statistics of the system components. 
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4.4 Density of quantum states 

The properties of a quantum system are governed by its Hamiltonian spectrum, whose form 
should also be important for quantum system equilibration. Because of the existence of 
the correspondence principle, one often compares quantum systems with classical ones. The 
equilibration of the latter depends on whether the system is integrable or not. Integrable 
classical systems do not equilibrate. For equilibration to happen, the system has to be 
nonintegrable. 

The notion of integrability is well defined for classical systems. If such a system, with 
n degrees of freedom, possesses n independent functions in involution (mutually Poisson 
commuting), then the system can be integrated up to quadratures [52]. This is known as 
Liouville theorem. 

Quantum integrable systems in n dimensions are often defined analogously by requiring 
the existence of n mutually commuting operators. However, there is no a theorem equiv- 
alent to the classical Liouville theorem, and there is nothing resembling the reduction to 
quadratures. In quantum mechanics, an integrable system is understood as that for which 
the spectral eigenproblem can be solved exactly [93]. The class of such exactly solvable 
quantum problems is very limited. 

Equilibration in closed classical systems is usually accompanied by the appearance of 
chaos [52]. The notion of quantum chaos is not well defined. One often resorts to the corre- 
spondence principle, telling that quantum chaos can exist, provided that the corresponding 
classical system is chaotic [94]. The latter requires that the system be nonintegrable. How- 
ever, classical chaos does not necessarily imply quantum chaos [95]. Quantum chaos is 
assumed to be related to the properties of the energy spectra. 

Gutzwiller [96] proposed that the spectra of integrable and nonintegrable quantum sys- 
tems should be qualitatively different, which would result in the qualitative difference of the 
density of states 

p{E) = Y,KE-E n ) . (58) 

n 

It has been suggested [96-100] that, generally, the spectra of quantum systems consist of two 
parts, regular and irregular, the regular part being almost equidistant, while the irregular 
part being random. Respectively, the density of states can be represented as a sum of two 
terms, regular and irregular [94,96,101,102], 

p(E) = p reg (E) + p irr (E) . (59) 

The regular term, in the semiclassical approximation, reads as 

Pre g (E) = J Q(E - U(v))[E - U {v)}^ dr , (60) 

where d is the space dimensionality and U(r) is an external potential. The semiclassical 
density of states ( 1601) is widely used for describing cold trapped atoms (see., e.g., [103]). If 
the external potential is a homogeneous function of degree a, such that 

U(Xr) = X a U(r) , 

then the regular term is 
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In particular, for a harmonic trap, when a = 2, one has 

Preg (E) (a = 2). 

The irregular term of the density of states can be represented as 

Pirr {E) Im ^ A n exp {z (s n - | /3 n ) } , (61) 

n 

where the sum is over semiclassical trajectories, S n is an action integral, and j3 n is the Maslov 
index [94,96,104]. The amplitude A n and the action integral S n are continuous functions of 
the energy levels. The irregular term of the density of states is a fastly oscillating function 
of energy. 



4.5 Statistics of level spacing 

The most evident difference between integrable and nonintegrable quantum systems is seen 
in the statistics of their level spacing. One considers the reduced spectrum spacing 

- m ' (62) 

where AE is the energy difference between the nearest energy levels, arranged so that AE 
be positive, hence s G [0, oo). 

For integrable quantum systems, the level spacing is random [99,105], being described by 
the Poisson distribution 

P(s) = e~ s . (63) 

The latter has the mean (s) = 1 and dispersion cr 2 = 1. The most probable energy gap is 
s P = 0. 

Nonintegrable quantum systems are characterized by the Wigner [106] distribution of 
level spacing 

P(s) = | s exp (- \ s 2 ) . (64) 
This has the mean (s) = 1 and dispersion <j 2 s = 0.273. The most probable energy gap is 

s w = = 0.798 . (65) 

The existence of the nonzero gap (|65|) is termed the level repulsion. 

Thus, integrable and nonintegrable quantum systems display essentially different statis- 
tics of their level spacing. Therefore equilibration can proceed differently for these systems. 
However, as is discussed in the previous subsections, integrable quantum systems can exhibit 
equilibration on average to a representative ensemble, though, may be, slower than nonin- 
tegrable systems. This is contrary to classical integrable systems that do not equilibrate at 
all. 



19 



5 Characteristic time scales 



5.1 Interaction time 

The process of equilibration of quantum systems from a strongly nonequilibrium state goes 
through several stages characterized by the corresponding time scales [7,62,82,107]. 

The shortest characteristic time is the time during which two particles interact with each 
other. For atomic gases, this time is of order tint ~ a s /v, with a s being a scattering length 
and v ~ l/ma s , atomic velocity, where m is mass. Thence, the interaction time reads as 

t int = ma 2 s . (66) 

For Hubbard-type models with the on-site interaction U, the interaction time is 

tint = JJ • (67) 

Generally, the interaction time is inverse to the typical strength of the strongest particle 
interactions. This time defines the first stage of relaxation 

< t < t int (interaction stage) , (68) 

when particles move being yet not correlated with each other. 

5.2 Local-equilibrium time 

After the interaction time, particles start developing mutual correlations that are growing till 
the local-equilibrium time that, for atomic gases, is of the order t\ oc ~ X/v, where A ~ 1/ ' pa 2 s 
is the mean free path and p is an average particle density. Then the local-equilibrium time 
is 

Tfi 

hoc = — ■ (69) 
pa s 

The interval of time between the interaction time and local-equilibrium time is the kinetic 

stage: 

tint < t < t[ oc (kinetic stage) , (70) 
when the motion can be described by kinetic equations. 

5.3 Heterophase time 

After the local-equilibrium time ti oc , the evolution enters hydrodynamic stage that can be 
subdivided into substages, when the system dynamics crosses a phase transition line. If the 
system starts from one phase, but is quenched into the conditions supporting another phase, 
then, after ti oc , there appear the nuclei of the new phase, which are also called clusters, or 
droplets, of this new phase, or mesoscopic topological defects. These are randomly intermixed 
in space with the regions of the former phase. Frenkel [108] called such systems heterophase. 
The nuclei of the competing phase can be of various shapes. In particular, they can be in 
the form of vortices [26,27,56]. Such a heterophase mixture exists during the heterophase 
stage 

tioc <t< t het (heterophase stage) , (71) 
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till the heterophase time thet, when the system becomes uniform, with the more stable phase 
filling the system volume. Equilibration passing through the heterophase stage has been 
observed, e.g., in experiments [20,22,24]. The heterophase stage can be of different duration, 
with thet being defined by the system properties [56,62]. In some cases, the heterophase time 
can be so long that the system can be treated as quasi-equilibrium during the whole stage 
( 17TT) . which happens for many condensed-mater systems [62,107,109-118]. 



5.4 Relaxation time 

The second substage of the hydrodynamic stage is that during which the system relaxes 
to an equilibrium or quasi-equilibrium state. This is the relaxation stage that lasts in the 
interval of time 

thet < t < t re i (relaxation stage) , (72) 

if there exists the heterophase stage (ITT!) . If the latter does not exist, then the relaxation 
stage starts from the local-equilibrium time (I691) . The heterophase stage may be absent, when 
the initial nonequilibrium conditions correspond to the same phase to which the system is 
assumed to relax, so that no transition line is crossed [119]. 

In numerical simulations for the Hubbard model, corresponding to atoms in optical lat- 
tices, the relaxation time was found to be 

trel = J ? (73) 

for both bosons [46,47] as well as fermions [48]. Introducing [18] the effective mass 

m * s • (74) 

where a is the distance between the nearest-neighbor sites, transforms the relaxation time 
f!73|) to the form 

t re i = Ira a 2 . (75) 

The effective mass, of course, can be very different from the atomic mass m. 

For quasi-one-dimensional Bose gases with local interactions, the relaxation rates were 
calculated by Mazets and Schmiedmayer [120]. The rate of populating the radially excited 
modes by pairwise atomic collisions, at T < u± was found to be 

r 2 = au ± exp ^- ^jr^j , (76) 

with the parameter 

a _ Na 2 s 
2lJ± 

Here iV is the number of atoms in the trap, l z = 1 / ^Jmuj z and ij_ = 1/ ^mu± are the effective 
trap lengths in the longitudinal and transverse directions, while u z and u± are the related 
trap frequencies. The relaxation rate, due to three-body collisions [120], is 

T 3 = 2a 2 co ± . (78) 
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The relaxation time is defined by the minimal inverse rate: 



t re i = min 




that is, by the maximal relaxation rate. Atomic correlations, described by the two-body and 
three-body correlation functions g n (0) at zero distance, can diminish [120] the relaxation 
rates to 

r n = r n # n (0) (n = 2,3), (80) 

resulting in the increase of the relaxation time. 

After the relaxation time, the system enters a quasi-equilibrium regime, lasting till the 
recurrence time. 



5.5 Recurrence time 

As is seen from Eq. fl35|) . the operator average is quasi-periodic, hence it always returns 
arbitrarily close to any given initial value after the recurrence time that can be estimated as 

2tt . . 

(81) 



AE 



where AE is the mean spacing between energy eigenvalues near E. The mean spacing is of 
the order AE ~ E/8e, with Se being the number of states in the interval AE. The entropy, 
as is known, is related to 5e as S ~ ln^g, thence 5e ~ exp(S'). For systems with N 3> 1, 
one has S ~ N and E ~ eN, with e being energy per atom. Therefore the recurrence time 
can be represented by the expression 

2vr Ar 

tree = Jn C (82) 

If the reduced energy e defines the interaction time t int ~ 1/e, then 

2ir N 

tree 6 tint ■ (S3) 

The recurrence time, for large N, can be extremely long. The quasistationary stage is 
associated with the time interval 

t re i < t < t rec (quasistationary stage) . (84) 

Generally, the typical relations between the characteristic times satisfy the inequalities 

tint < t[ oc < thet < t re [ < t rec . (85) 

Numerical values of these time scales, for condensed-matter systems and atomic gases, are 
given in Refs. [7,62]. 
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6 Nonisolated quantum systems 
6.1 Quantum system in a bath 

If a finite quantum system is not isolated, one needs to treat the coupled pair of the system 
and its surrounding, called bath. The composite system is characterized by the space of 
microstates 

n AB = n A 0n B , (86) 

being the tensor product of the Hilbert space 7i A , related to the studied system, and of the 
Hilbert space Hb, related to the bath. The dimension of the composite space is 

dimH AB = d A ds , (87) 

where the factor space dimensions are 

d A = dimH A , d B = dimH B ■ (88) 

The total Hamiltonian is the sum 

H AB = H A (g) l B + l A (g) H B + H int , (89) 

where 1a and 1^ are the unity operators defined on % A and Hb, respectively. To simplify 
the notation, one usually omits the unity operators, writing 

H AB = H A + H B + H int , (90) 

which implies Eq. ( 1891 . 

The statistical operator of the composite system is denoted by p A B (t) ; which can be pure 
or not. One introduces partial statistical operators 

p A (t) = Tr B p AB (t) , p B (t) = Tr A p AB (t) , (91) 

in which Tr A means the trace over the space H, A and Ttb, the trace over ~H.b- 

Observable quantities are related to the system defined on % A . The expectation value of 
an operator A of a local observable, acting on T-i A) is 

(A(t)) = Tr AB p AB (t)A = Tr A p A (t)A , (92) 

where Tr AB is the trace over the space H A b- 



6.2 Distance between statistical operators 

If one is interested in the equilibration on average, one needs to consider the time-averaged 
statistical operator of the composite object, 

p AB = lim - / p AB {t) dt , (93) 

T->00 r J Q 



23 



and the time-averaged partial statistical operators 

p A = lim - f p A (t) dt 



p B = lim - / p B (t) dt . (94) 



T— >oo T 







The difference between two statistical states, described by the operators p\ and p 2 , defined 
on the same Hilbert space H, can be characterized by the Hilbert-Schmidt trace distance 



dist [pi, p 2 ] = Tr^y (pi - p 2 ) 2 . (95) 
Respectively, the time-averaged distance is 

dist[p!, p 2 ] = lim - / dist[pi(t), p 2 (t)] dt . (96) 

It has been shown [121-123] that 



dist [p A (t), p A ] < yJd A (TX?B) ■ (97) 

In the proof, one uses the Jensen inequality, telling that for a concave function, for which 
d 2 f(x)/dx 2 < 0, one has f(x) < f(x), and the inequality 



Tr^y (pi - p 2 ) < y d A Tr (pi - p 2 ) . 
If the bath is large, with a great number of uniformly distributed states, such that 

Trp| « i- , (98) 

then the system equilibrates on average in the sense of the small trace distance (1971) . 

6.3 Decoherence induced by environment 

For the problem of equilibration, a large bath, representing an environment in thermody- 
namic limit, is of special interest. The standard definition of thermodynamic limit for a 
system of iV particles in a volume V reads as 

N 

N — t- oo , V — > oo , — — > const . 

This definition, however, is not applicable for nonuniform quantum systems confined in an 
external potential, such as trapped atoms, since the system volume V is not well defined in 
this case. A general definition of thermodynamic limit, which is valid for arbitrary systems, 
including confined ones, is given as follows [18,103]. For any extensive observable quantity 
An =< A >, corresponding to a system of N particles, it should be: 

N — > oo , A N oo , — > const . 
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As is obvious, for a uniform system, where the volume V is well defined, the latter definition 
reduces to the standard one. 

But we are interested in thermodynamic limit for a bath of volume V B and described 
by an effective number of degrees of freedom Nb, while the considered quantum system is 
kept finite, with a finite number of atoms N. Then thermodynamic limit for the bath can 
be defined [124] in the standard way 

N R 

N B -» oo , V B ->• oo , — > const . (99) 

Vb 

The role of the bath can be played by other parts of the same physical system, where a finite 
part of it is fixed. In that case, thermodynamic limit is applied to the whole system [125]. 
The effect of decoherence assumes that at long times the average (192]) yields 

lim lim (A{t)) = Y J PnnKn, (100) 

n 

where the bath thermodynamic limit ( 19 9 p is implied. The thermodynamic and time limits 
do not commute with each other. Decoherence, as is defined in Eq. ( llOOp . also means 
equilibration, according to definition (J1J. 

Decoherence, as a general phenomenon, is supposed to hold for any operator A from the 
algebra of local observables. Condition ( 1100}) may be trivially valid for some operators that 
commute with the system Hamiltonian Ha, which has nothing to do with decoherence. 

The reduction of an operator average to the diagonal sum fllOOj) is called decoherence or 
dephasing [126,127] because of the evident reason. In quantum mechanics, the nondiagonal 
terms characterize interference effects that are typical of coherent processes. In radiating 
systems, such nondiagonal terms are responsible for the appearance of coherent radiation 
[71-80]. Hence, the disappearance of the coherent terms is nothing but decoherence. 

The properties of the environment are often characterized by random distributions. For 
examples, expanding the wave function of the composite system (system plus bath) over 
a basis, with the real and imaginary parts of the expansion coefficients being independent 
Gaussian random variables with zero mean, one comes to the reduced statistical operator 
for the system, having canonical form [128,129]. 

Assuming a random distribution of the wave-function expansion coefficients is equivalent 
to assuming infinite number of degrees of freedom of environment. It is therefore possible, to 
impose the condition of randomness on the wave function of the considered system, without 
explicitly invoking a bath [130-132]. 

Similarly, equilibration and decoherence occur, if one assumes that the system correlation 
functions exponentially fastly decay [133,134]. This assumption is analogous to the random- 
phase averaging and can be interpreted as a random influence of environment. 

7 Quasi-isolated quantum systems 
7.1 Absence of complete isolation 

The assumption of the possibility of an absolute system isolation is, evidently, an idealization 
of the real situation. In reality, no one system can be completely isolated. There always 



25 



exists some random influence of environment, though may be weak. Moreover, the notion 
of absolute isolation is self-contradictory, since in order to state that the given system has 
been isolated during a period of time, it is necessary to accomplish a series of measurements 
with this system for proving that the system has really been isolated during that time 
span [60,135,136]. But each measurement, even if it is a nondemolition measurement [137], 
disturbs the system, making it not isolated. Measuring devices play the role of environment 
acting on the system. And, vice versa, the influence of environment is analogous to the 
action of measuring procedures. 

In computer modeling, the role of external noise is played by numerical errors accompany- 
ing calculations. With time, there is the accumulation of errors, analogous to the increasing 
influence on the system of random environment. Thus, a complete isolation is merely a 
theoretical abstraction not existing in nature. There can exist only quasi-isolated systems. 

7.2 Properties of quasi-isolated systems 

To formulate rigorously the notion of quasi-isolated systems, let us proceed to mathematical 
definitions. Let Ha be a system Hamiltonian defined on the Hilbert space Ha and let Hb, 
defined on Hb, be the Hamiltonian of surrounding. The total space of microstates for the 
composite object "system plus environment" is the tensor product 

n A B = nA0n B ■ (101) 

The Hamiltonian of the composite object is 

Hab = Ha + H B + H int , (102) 

where Hi nt characterizes interactions between the system and environment and where the 
unity operators, for the brevity of notation, are omitted. 

It is reasonable to call a system quasi-isolated, when the matrix representation for the 
operators of observable quantities, in the basis of the eigenvectors of the system Hamiltonian 
Ha, is not changed under the influence of the environment, so that this matrix representation 
coincides with the matrix representation of these operators in the basis of the eigenvectors 
of the total Hamiltonian (11021) . This means that Ha is to be the integral of motion with 
respect to Hab, 

[H a , H ab ] = . (103) 

Condition (11031) is the mathematical definition of a quasi-isolated system. The latter also 
implies the commutation relation 

[H A , H int \ = . (104) 

The property of quasi-isolatedness defines nondestructive, or nondestroying, action of en- 
vironment on the system, similarly to the minimally disturbing action of measurements 
[60, 135,136]. This property should not be confused with the definition of nondemolition 
measurements for a given observable [137], where the commutation of the operator of the 
measured quantity with the total Hamiltonian is required: [A, Hab] = 0. 
When the system Hamiltonian satisfies the eigenproblem 

H A \n) = E n \n) , (105) 
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and the system is quasi-isolated, then the eigenproblem for Hab has the form 

H AB \nk) = (E n + e nk )\nk) , (106) 

where the vector \nk >= |n > (g>\k > pertains to the space %ab- 

For the statistical ensemble {Hab, PAB(t)}, the evolution of the statistical operator is 
given by the law 

PAB(t) = U A B(t)p A B(0)Ui B (t) . (107) 

If Hab does not explicitly depend on time, then 

U AB (t) = exp(-iH AB t) . (108) 
As an initial condition for the evolution of the statistical operator, we accept 

Pab{0) =Pa{0)6§Pb{0) > (109) 



which means that the interaction between the considered system and the bath is switched on 
at if: = +0. This initial condition is not principal, but it is assumed here merely for simplicity. 

7.3 Equilibration of quasi-isolated systems 

We are interested in the temporal behavior of operators of local observables, associated with 
the system. Say, we consider the average of an operator A defined on %a- For this average, 
we have the expression 

(A{t)) = Tt AB p AB {t)A = Tr A pA{t)A , (110) 

in which 

p A (t) = Tr B p AB (t) . (Ill) 
In the matrix representation, this yields 

(A(t)) = J2pL(t)A nm , (112) 

ran 

where 

pi n {t) = (m\p A (t)\n) = J2( mk \pAB(t)\nk) . (113) 

k 

In view of the evolution operator f)108p . Eq. (I113P can be written as 

Pn\n(t) = P^nk exp{-i(u mn + E mnk )t} , (114) 
k 

with the notation 

pi* k = (mk\p AB (0)\nk) (115) 

and 

^mn = E m E n , £mnk = &m,k &nk ■ (H6) 
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Thus for the average fill 01) . we get 

(Mt)) = Yl Yl P™ B nk A nm exp{-i(u mn + S mnk )t} . (117) 
mn k 

Taking into account the initial condition (I109p yields 

Pi* k = (m\p A (0)\n) (k\p B (0)\k) . (118) 
The latter, employing the notation 

Pmn = (m\p A (0)\n) , p k = (k\p B (0)\k) , (119) 

reduces to 

Pink = PmnPk ■ (120) 

The forms p mn and p k satisfy the normalization conditions 

J2pnn = l, J>fc=l- (121) 

n k 

Introducing the decoherence factor 

D mn (t) = J2Pke- i£mnkt , (122) 

k 

we come to the expression 

Pmn(t)A nm D m n(t) j (123) 

mn 

in which 

p mn (t) = pmne-™™ 1 . (124) 
Separating in summation (11231) the diagonal and nondiagonal terms gives the average 

(A(t)) = ^ ^ PnnA nn + ^ p mn (t)A nm D mn (t) , (125) 
n m^n 

where we take into account that, according to definition (jllfip . 

Sunk = . (126) 

Introducing the density of states 

k 

satisfying the normalization condition 

Pmn (e) de = 1 , (128) 

oo 
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transforms the decoherence factor (I122p to 



Dmnit) = / Pmn (e)e- i£t de . (129) 

J — oo 

For the diagonal elements, we have 

Pnn{e) = ^2p k 5(e) , D nn (t) = 1 . 
k 

Therefore the decoherence factor is nontrivial only for the nondiagonal elements of Eq. (11291) . 

Assume that the surrounding environment is large, such that the spectrum of Hab, with 
respect to the index k, can be treated as continuous and the summation over k can be 
replaced by integration. Then the density of states (11271) may be considered as a smooth 
function of e. In that case, we can invoke the Riemann-Lebesgue lemma [138] telling that, if 
p mn (e) is a measurable L 1 -integrable function, then the time limit of Eq. (I129p . for m ^ n, 
is zero: 

lim D mn {t) = (m ^ n) . (130) 

t— >OD 

Hence for average (11251) . we obtain 

Mm(A(t)) = Vp„nAm • (131) 
t— >oo ' 4 

n 

Therefore, a quasi-isolated system in a large environment equilibrates. This has been 
proven here for an arbitrary quasi-isolated system, which makes the given proof principally 
different from the considerations, where the property of quasi-isolatedness was not used, but 
instead, concrete systems were studied. For instance, equilibration for quadratic Hamiltoni- 
ans was demonstrated [124,134]. Also, we have considered equilibration in the strict sense 
of definition (j3J), but not equilibration on average [86-92], as described in Sees. 4.2 and 4.3. 

7.4 Examples of decoherence factors 

To illustrate explicitly admissible forms of the decoherence factor (11291) . let us consider 
several concrete cases. 

When the density of states enjoys the Gaussian form 

Pmn(e) = — = exp (- ' ( 132 ) 

then the decoherence factor is also Gaussian, 



D mn (t) = exp I - % t 2 ) . (133) 



For a uniform density of states 

= 2 T 



p mn {e) = 7— 6(7 mn - e)Q(j mn + e) , (134) 
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where 9(e) is a unit-step function, the decoherence factor is 

n(7 m r. 

Ifmnt 



A»W = 5^ . (135) 



Assuming the Poisson distribution 



27 m n V 7- 
yields the Lorentz-type decoherence factor 



p mn {e) = ^— exp ( - ) , (136) 



DUt) = TT(wP ' < 137) 

Conversely, for the Lorentz density of states 

p ™" (£) = ' (138) 

the decoherence factor is exponential, 

D mn (t) = e-^»* . (139) 

In the above examples, the density of states is represented by monotonic nonincreasing 
functions of e. For nonmonotonic functions, the situation is similar, just formulas can become 
a bit more cumbersome. For instance, if we take the density of states in the Wigner form 

7r|£| / HE 1 



Pmn{e) = — expf-^— ) , (140) 

/ mn \ I ran , 



then the decoherence factor reads as 



D mn (t) = 2 J™ xe~ x2 cos x) dx . (141) 



The integral 

2 J xe~ x2 cos(qx) dx = 1 — qF (J^j 
is expressed through the Dawson function 

F(z) = e- z2 f e x2 dx . 
Jo 

And for the decoherence factor we find 



n (/) i - ^p tF (ln^.) , (ij2) 

V 71 " V v 71 ' 



At large variable, the Dawson function behaves as 

F(z)*±- + ± (*->oo). 
Consequently, the decoherence factor decays as 

D mn (t) ~ - — (t^oo). (143) 

Iran 

It is seen from these examples that the typical decoherence time is 

Uec = ~ , (144) 

7 



where 7 is an average r y mn . 



30 



7.5 Occurrence of rare events 



The property of asymptotic decoherence (11301) and, consequently, of the equilibration limit 
(I13ip . are based on the use of the Riemann-Lebesgue lemma assuming that the density 
of states p mn (e) is a measurable L 1 integrable function. By its definition, the density of 
states is a non-negative function that is normalized as in Eq. (I128p . Hence, it is always L 1 
integrable. But it may happen that it can contain nonmeasurable parts. Then the properties 
of decoherence and equilibration are not guaranteed. 

For instance, the density of states can be formed of two parts, 

Pmn{e) = XiP^ie) + AaiC(e) , (145) 

in which the first term is a measurable function leading to a relaxation process accompa- 
nied by decoherence and equilibration, while the second term is a nonmeasurable function 
resulting in decoherence oscillations, and where 

Ai + A 2 = 1 (Ai > , A 2 > 0) . (146) 

Thence the decoherence factor is also a sum 

D mn {t) = \iD&{t) + X 2 DZXt) (14T) 

of the related terms, the first of which satisfies property (11301) . 
Let the second term of the density of states have the form 

= \ ~ «i) + 6 ( £ + > ( 148 ) 

3 

where atj are real quantities. Then the oscillating part of the decoherence factor is 

3 

Therefore the level of decoherence will oscillate, being a quasi-periodic function characterized 
by the periods 

t s = — . (150) 

atj 

Respectively, the operator average (11231) does not decohere and the system does not equili- 
brate. However it may look as equilibrium on average, keeping in mind temporal averaging. 

If the parameters atj are small, nonequilibrium oscillations will be rare. But if they are 
strong, they cannot be neglected, since their presence can essentially influence the system 
properties. Such rare, but strong, events can support the existence of heterophase fluctua- 
tions [62,107,109-118]. 

8 Equilibration under nondestructive measurements 
8.1 Definition of nondestructive measurements 

As has been mentioned above, measurements accomplished with the considered system act 
similarly to the action of environment. The main difference is that the measurement proce- 
dures, generally, depend on time, which makes the consideration more complicated. 



31 



Let the system be characterized by a Hamiltonian Ha, defined on a Hilbert space Ha, 
and a measuring device be described by a Hamiltonian Hd, defined on a Hilbert space Tip. 
These Hamiltonians Ha and Hd are assumed to be time independent. 

The total space of microstates is 

U = n A (g)H D . (151) 
The total Hamiltonian can be written as 

H = H A + H M , (152) 

where the measurement Hamiltonian 

Hm = Hd + H int (153) 

is the sum of the Hamiltonian Hp of the measuring device and of H int characterizing the 
interactions between the studied system and the measuring device. The system Hamiltonian 
Ha and the device Hamiltonian Hd do not explicitly depend on time, but the interaction 
Hamiltonian Hi nt is a function of time. Consequently, the measurement Hamiltonian is 
an explicit function of time, H M = H^(t). As usual, we omit in the formulas the unity 
operators, keeping in mind that the rigorous expression for Hamiltonian f)152p would be 

H = H A (g) Id + U (g) H D + H int . (154) 

We consider measurements that minimally disturb the system [60,135,136], because of 
which they are named nondestructive. The observed system has to be quasi-isolated from 
the measuring device. A measurement is called nondestructive, if the system Hamiltonian, 
as well as the device Hamiltonian do not disturb the properties of each other, being the 
integrals of motion: 

[H A , H]=0 , [H D , H] = . (155) 

It follows from Eq. (11551) that 

[H A , H m ] = , [H D , H m ] = . (156) 

If the eigenproblem for Ha has the form 

H A \n) = E n \n) , (157) 

then the measurement Hamiltonian satisfies the eigenproblem 

H M (t)\nk) = e nk (t)\nk) , (158) 

where \nk >= \n > ®\k >, with {\n >} being a basis in %a and {\k >}, a time-independent 
basis mH D - 

The statistical operator for the composite system is 

pit) = U(t)p(0)U + (t) , (159) 

in which the initial condition 

m=PA(0)6dp D (0) (160) 
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is assumed. This means that measurements start at the moment of time t = +0. 

In view of eigenproblem (I158p , it follows that the measurement Hamiltonian satisfies the 
Lyappo-Danilevsky condition 



H M (t), f H M (t')dt' 
Jo 



. 



(161) 



Since the system Hamiltonian Ha does not explicitly depend on time, the previous property 
f )16ip gives 



H(t), [ H(t')dt' 
Jo 

Therefore the evolution operator has the form 



. 



U[t) = exp l-ij H(t') dt' 
With the total Hamiltonian (I152p . this becomes 

U(t) = exp l-iH A t - i [ H M (t') dt' 



(162) 



(163) 



(164) 



(165) 



(166) 



8.2 Temporal dependence of observables 

Observable quantities are represented by operator averages. Let us consider the evolution of 
an observable quantity represented by an operator A acting on T-La- This observable is given 
by the average 

(Alt)) = Tr w p(t)i = Pmnk{t)A nm , (167) 



In view of eigenproblems f ll57p and fll58p . we have 

H\nk) = (E n + e nk (t))\nk) . 
Therefore the action of the evolution operator on the vector \nk > yields 

U{t)\nk) = exp \-iE n t - i [ e nk {t') dt'X \nk) . 



rank 



in which 



Pmnk{t) = (mk\p(t)\nk) . 
Under the evolution operator (I164p . we have 

-t 



Pmnk(t) = Pmnk(0) exp j-iu; mn t - % J E mnk {t') dt' 



with 



Pmnfe(O) = (mk\p(0)\nk) 
and where the notation is used for the energy differences 



£mnk(t) = £mk(t) ~ ^nk(t) W ' 



= E - E 

mn — J - / m - t -'n 



(168) 

(169) 
(170) 
(171) 
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Notice that 

Snnk(t) = . (172) 

Because of the initial condition ( I160p . Eq. (1 170 ft reduces to 

Pmnfc(O) = pmn{0)Pk , (173) 

where 

Pmn (0) = (m\p A (0)\n) , Pk = (k\p D (0)\k) . (174) 
The normalization conditions 

J^Pnn(0) = l, ^ Pk = l 

n k 

are valid. In that way, Eq. (ll69p becomes 

Pmnk{t) = pmn{t)PkeXp \ -% / S mnk {t') dt' \ , (175) 



with 

Pmnit) = p mn (0) exp(-iu mn t) . (176) 
Introducing the decoherence factor 

D mn (t) = ^p fc exp j-^y £mnk(t') dt'^j (177) 
allows us to represent average f ]16Tf) in the form 

Pranif) A nm D mn (t) . (178) 

mn 

Taking into account properties (11711) and (11721) . we can rewrite the above average as 

(A(t)) = ^PnnAnn + Pmn(t) A nm D mn (t) , (179) 

with p nn = pnn(0). This expression is similar to the average for a quasi-isolated system, but 
with a different decoherence factor. 

8.3 Specification of measurement procedure 

To better understand the structure of the decoherence factor (j!77p . we need to specify the 
action of the measurement Hamiltonian. Let this Hamiltonian, describing M measurements, 
have the form 

M 

H M (t) = J2Xjfj(t), (180) 

3=1 

where fj(t) is a real function and Xj is a self-adjoint operator, not explicitly depending on 
time. This Hamiltonian satisfies the Lyappo-Danilevsky condition (11611) . According to Eq. 
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, the operator H M and, hence, Xj, possess the eigenvectors \nk >. That is, we can 
write 

Xj\nk) = Xj n k\nk) , (181) 

where Xj nk is a real number, characterizing the measurement impact. As an explicit illus- 
tration, we can define Xj as a spectral resolution over projector operators: 



Xj = ^x jnk P nk (P nk = \nk)(nk\j . 



nk 

The eigenvalues of Hm acquire the form 

M 



e n k{t) = J2 x ^kfj{t). (182) 

3=1 

Respectively, for the transition frequency, defined in Eq. f)17ip . we get 

M 

£mnk{t) = Aj mnk fj(t) , (183) 

3=1 

where 

Aj mnk = Xjmk Xjnk ■ (184) 

Thence the decoherence factor (11771) becomes 

f M 1 

D mn (t) = ^Pk exp < -i ^ Ajmnk¥j(t) y , (185) 

k { j=l J 



with the phase function 



<Pj(t) = / /i(0 dt' ■ (186) 



Being interested in the long-time behavior of the decoherence factor, let us assume that 
the phase function, after the last measurement, occurring at the time tu-, does not depend 
on the measurement index j, 

ifj{t) = <p(t) (t > t M ) . (187) 
Then, after t > t^j, the decoherence factor takes the form 

D mn (t) = 22p k exp{-iA mnk Mip(t)} , (188) 

k 

in which the definition of the mean impact 

- 1 M 

A mnk = ^ ^ Aj mnk (189) 

3=1 
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is used. Finally, introducing the distribution of the measurement impacts 

Pmn{x) = ^PkS (x - A mnk ) , (190) 

k 

we come to the decoherence factor 

D mn (t) = J p mn (x) exp {-ixM<f(t)} dx . (191) 

The measuring device is a macroscopic object, because of which the spectrum of the 
measurement Hamiltonian, with respect to the index k should be treated as continuous. 
This means that the notation has to be understood as integration over k. Consequently, 
the impact distribution (11901) can be considered as a measurable function. By its definition, 
it is also L 1 integrable. Therefore, by the Riemann-Lebesgue lemma, 

Dmn(t) -> (Mip(t) -> oo) . (192) 

That is, nondestructive measurements lead, at large times, to the system decoherence, pro- 
vided that either the number of measurements M or the phase function ip(t) increases with 
time, so that their product also tends to infinity. 



8.4 Decoherence under nondestructive measurements 

To obtain an explicit form of the decoherence factor, we have to specify the impact distri- 
bution (I190p . Taking for the latter a Gaussian distribution 

p mn (x) = —=t exp (- \ (193) 

V^TTOmn I ACT mn ) 

yields the decoherence factor 

'2 



D mn (t) = exp <j - MVW } ■ (19D 
In the case of the Lorentz distribution 



Or, 

n (x 2 + a 2 mn ) 



(195) 



we have the decoherence factor 

D mn (t) = exp{-a mn M<p(t)} . (196) 

Let us consider two opposite types of measurements, instantaneous and continuous. For 
instantaneous measurement, we have 

= *(* - *i) , i Pj (t) = e(t-t j ). (197) 

Hence, 

l p(t) = cp j (t)\ t>tM = l. (198) 
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And continuous measurement implies 



fj{t) = 1 , v^t) = t . (199) 

In any case, decoherence happens when the number of measurements M increases [60,135,136]. 

If measurements are accomplished during all period of observation, the number of mea- 
surements is proportional to the time elapsed. This can be formalized by the relation 

t 1 - 

To summarize the results, we omit, for brevity, the indices m and n, marking instead 
the form of the impact distribution p(x), Gaussian (11931) or Lorentzian (1195ft . and the type 
of the measurement, instantaneous or continuous. We assume relation ( 1200j) and define the 

decoherence time 

t dec EE J — . (201) 

V o 

Then we obtain the decoherence factors, under instantaneous measurements, for the 
Gaussian distribution, 



1 / t 



2 




Z^(t) = exp } } . (202) 

and for the Lorentzian distribution, 

DLt(t) = exp (- . (203) 

While in the case of continuous measurements, we find for the Gaussian distribution, 

(204) 

and for the Lorentzian distribution, 

^L t (*) = expj- (^y| . (205) 

Continuous measurements result in much faster decoherence than instantaneous ones. 

When decoherence leads to the disappearance of the second term in the operator average 
( I179p . this also implies equilibration. Therefore, equilibration of a finite quantum system can 
be achieved by the measurement procedure. As has been emphasized above, the notion of 
complete isolation is self-contradictory, since to state that the system is isolated requires ac- 
complishing measurements confirming the fact of isolation [60,135,136]. But measurements, 
even the nondestructive ones, produce decoherence in the observable quantities. In the long 
run, decoherence, induced by measurements, results in the system equilibration. 
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8.5 Existence of time arrow 



The real time in the world, as is known, flows in one direction, from the past to future. 
The real time arrow cannot be reversed. But the microscopic equations of motion are 
time-symmetric, allowing one to reverse the direction of time. The problem why this time 
symmetry becomes broken, when passing from the microscopic equations to the observable 
quantities, is the long standing puzzle attracting the attention of many. It is not our aim to 
give here a discussion of various approaches to this problem. But we want to stress that this 
problem finds a simple resolution in the frame of the equilibration theory of quasi-isolated 
systems. 

The explanation for the existence of the unidirected time arrow lies in the fact that 
completely isolated systems do not exist in nature [139-141]. In the best case, they can 
be quasi-isolated, interacting with environment that, even if it does not destroy and does 
not essentially disturb the system, anyway, influences the temporal behavior of observable 
quantities. Even if the influence of environment is so weak that it can be neglected, the 
considered system is subject to measurement procedures. An absolutely isolated system is 
an abstraction that cannot be observed. As soon as one passes from time-reversible micro- 
scopic equations of motion to the consideration of observables, on has to take into account 
that the studied system can only be quasi-isolated. And the statement that the system is 
quasi-isolated necessarily presupposes that this fact is to be confirmed by measurements. In 
any case, whether there exists environment acting on a quasi-isolated system or the latter is 
subject to at least nondestructive measurements, there appears decoherence and equilibra- 
tion. Each of the latter phenomena influences the time dependence of observable quantities, 
making their evolution irreversible. 

The existence for any quasi-isolated system of attenuation, caused either by environment 
or by measurements, or by both, leads to equilibration in the strict sense of limit but 
not merely in the sense of time averaged quantities. Even infinitesimally small such an 
attenuation principally changes the situation, resulting in the existence of time arrow. 

In order to emphasize the importance of even an infinitesimally weak nonisolatedness, 
let us consider two limits: One limit is when the measurement dispersion a tends to zero, 
which corresponds to neglecting the influence of measurements or environment. In the case 
of the environment influence, the role of a is played by the environment attenuation rate 7. 
Another limit is when time tends to infinity. For the decoherence factor, in the first case, 
we have 



HmD(t) = 1 . 



(206) 



While in the second case, 

lim D{t) = . 



t— >oo 



(207) 



These properties lead to the noncommutativity of the limits, as far as 



a— >0 t—>oo 



lim lim D(t) = , 




while reversing the limits, we get 



t—>-00 CT-S-0 



lim lim D(t) = 1 . 



(209) 
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This noncommutativity of the limits for the decoherence factor results in the noncom- 
mutativity of the limits for the observable quantities: 



lim lim (A(t)) ^ lim \im(A(t)) . (210) 

CT— >(} t— >00 t— >00 CT— >0 

The noncommutativity of these limits is a pivotal property of quasi-isolated systems, ex- 
hibiting the existence of time arrow [139-141]. 

We may also notice that, if we would like to formally consider negative time, then the 
exponentials in the decoherence factors f ll3QI> and ( 12031) should be treated as depending on 
\t\. Therefore, in all cases the average < A(t) > equals < A(—t) >. That is, formally 
the observable quantities are symmetric with respect to time inversion, in agreement with 
microscopic equations of motion. However, any nonequilibrium statistical state of a quasi- 
isolated system always tends with time to an equilibrium or, at least, to a quasi-equilibrium 
state. This tendency explains the existence of time arrow. 



9 Summary 

Finite quantum systems, because of their important and widespread technological applica- 
tions, have become the topic of intense investigations, both experimental and theoretical. 
Nonequilibrium properties of such systems are of special interest, being less understood than 
equilibrium ones. The principal question is whether finite quantum systems equilibrate and 
if so, when and how? 

Cold trapped atoms provide a very convenient laboratory for the study of nonequilib- 
rium phenomena. The most interesting experiments with nonequilibrium trapped atoms are 
reviewed in the present survey, and the related numerical simulations are discussed. 

Finite quantum systems, strictly speaking, cannot be absolutely equilibrium, as far as 
their statistical states are quasi-periodic. But, the studies show that isolated integrable 
as well as nonintegrable quantum systems can equilibrate on average to a quasi-equilibrium 
state with a rather long lifetime, since the recurrence time is usually very long. Equilibration 
to such quasi-stationary states can be faster or longer depending on the system integrabil- 
ity The resulting quasi-stationary states are, generally, described by representative Gibbs 
ensembles. In particular cases, these can be microcanonical ensembles. 

Nonisolated quantum systems, interacting with a macroscopic bath, equilibrate in the 
strict sense. Equilibration is accompanied by decoherence. The time dependence of the 
decoherence factor is defined by the density of states that, in turn, depends on the system 
parameters. The relaxation can be exponential or Gaussian [142]. 

Actually, completely isolated finite systems do not exist, but there are only quasi-isolated 
systems. Any finite system is always influenced by its environment, maybe weakly, but, 
anyway, noticeably. Moreover, the concept of isolated systems is self-contradictory, since in 
order to state that a system is isolated during a period of time, it is necessary to prove this 
by a series of measurements, whose influence makes the system not isolated. Such quasi- 
isolated systems can equilibrate in the strict sense. In some cases, there can exist rare events 
making the resulting state quasi-equilibrium. The latter is equilibrium on average. 

In any case, quasi-isolated systems do relax to an either equilibrium or a quasi-equilibrium 
state. This irreversible behavior explains the existence of time arrow. 
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In this review, the equilibration processes have been described, starting from a strongly 
nonequilibrium state and tending to a steady state. Another problem of great interest is 
the opposite process of the development of nonequilibrium states from the given equilibrium 
one. For example, acting on trapped Bose atoms by external alternating fields [18,143,144], 
it is possible to generate dynamic transitions from the ground-state Bose-Einstein conden- 
sate to a vortex superfluid and then to a turbulent superfluid [145-147]. Investigating such 
transformations from a given equilibrium state to a highly nonequilibrium one would allow 
for the better understanding of the development and properties of nonequilibrium states in 
finite quantum systems. 

Nonequilibrium effects in finite quantum systems can find numerous applications, ranging 
from various electronic devices to quantum computers. As has been noticed long tome ago 
by Schrodinger [148] and Bohr [149], even the functioning of humans may require, for their 
correct description, the necessity of employing the theory of finite quantum systems. For 
instance, this could be necessary for understanding the functioning of genes and brain. In 
turn, understanding the brain activity would allow for the creation of thinking quantum 
systems [150]. 
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